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Abstract 

It is shown that the scattering of electromagnetic waves by Langmuir ones taking into the account 
the electric drift motion of particles, is the most intense nonlinear process and should be responsible 
for the formation of radiation spectrum of radio pulsars. Performed analysis indicates that in the 
case of existence of inertial interval formation of stationary spectra is possible. Analysis of linear 
mechanisms of wave generation and silk allow to conclude that only possible stationary solution 
has spectral index a = —1.5. Obtained spectrum is in good agreement with observational data. 
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I. INTRODUCTION 



Aim of the present work is to study nonlinear processes responsible for the formation 
of the radio emission spectrum in the pulsar magnetosphere. The spectra of pulsar radio 




0. 



Above 100 MHz observed 



emission has been intensively studied by observers [1, 
intensity / of the majority of pulsars can be described by a simple power law / ~ v a with 
average value of spectral index (a) = —1.8 ± 0.2 2]. Spectral index of individual pulsars 
vary from a = +1.4 to a = —3.8. Some pulsars has the spectra which can be modelled by 
the two power law Q,!^, but it seems that single power law is a rule and the two power law 
spectrum is a rather rare exception Q]. 

From our point of view pulsar radio emission is generated and formed in the magneto- 
sphere. Therefore properties of the radiation is mainly determined by the physical processes 
in the magnetosphere plasma. This approach allows to explain other main characteristics of 
pulsar radiation, such as polarization properties jfj, mailings L 7|, micro impulses mode 
switching |2| and so on. 

In the presented paper we study nonlinear processes in the magnetosphere plasma that 
should be responsible for the formation of the pulsar radio emission spectra. Mathematical 

nfi 

methods used for analysis is based on the weak turbulence theory (see, e.g., 10. 111). 

The paper is organized as follows: the model of the magnetosphere as well as linear 
eigenmodes of e~e + plasma and linear mechanisms of their generation is discussed in Sec. 
2. Nonlinear processes are studied in Sec. 3. In Sec 4 stationary solutions fot the spectra 
are obtained and analyzed. 



II. LINEAR MODES AND THEIR GENERATION IN PULSAR MAGNETO- 
SPHERE 



According to the standard model of the pulsar magnetosphere [12J, [13J] relativistic e~e + 
plasma that is pierced by electron beam moves along open field lines of pulsar magnetic field 
B . Magnetic field of pulsar is supposed to be dipole. Under open field lines we mean the 
field lines that crosses the surface of light cylinder - the surface where the velocity of field 
line solid rotation reaches the speed of light c. It is supposed that electric field directed along 
the open field lines is generated at the surface of pulsar (l^ . that pulls out electrons of the 
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primary beam from the surface of the star. The particles moving along magnetic field lines 
generate 7 rays, which on its turn if the energy e 7 > 2m e c 2 generate electron-positron pairs 



13j. This process is stopped when the initial electric field becomes screened. The particle 



distribution function that has the form presented in Fig. ^ can be presented as: 

f = f P + ft + f b , (1) 

where f p describes the bulk of the plasma, ft is the distribution function of elongated to the 
plasma motion tail particles and /& corresponds to the particles of the primary beam. The 
distribution function (P) is one dimensional. In Fig. ^ solid and dashed lines correspond to 
electrons and positrons respectively. The shift of the distribution functions is caused by the 
existence of the primary beam. It is supposed that 

rip'jp « n^t « ^rib'jb, (2) 

where 7 P ,7t, 76 are Lorentz factors and n p ,n t , rib concentrations of bulk, tail and beam 
particles respectively. For typical pulsars jb ~ 10 6 — 10 7 , n° ~ 10 11 , j t ~ 10 4 , n° t ~ 
10 13 — 10 14 , 7 P ~ 10 rip ~ 10 17 . Index indicates that the values belong to the surface of 
the pulsar. 

The condition of quasi neutrality yields 

A 7 = 7+ - 7- = J Uld 3 p - J /_ 7 d 3 p, (3) 

here indexes ± are related to positrons and electrons respectively. f± are normalized such 
that / f±d 3 p = 1. 

In spite of the fact that A7 <C 7-t this small difference plays important role in explanation 
of polarization properties of pulsar radiation. The radiation of the most pulsars has consid- 
erable part of circularly polarized radiation. It was shown |(| that circular polarization can 
occur only in small angle 9 between wave veqtor k and B : 

9 2 < — A 7 « 1. (4) 

UJ B 

Consequently, the observer receives the radiation that is formed in small angle with respect 
to Bo- So, for pulsar radio emission the waves propagating nearly along the magnetic field 
are important. This agrees with linear theory of e~e + plasma in strong magnetic field: 
analysis provides that the wave generation takes place only almost along the magnetic field 
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Due to the absence of gyrotropy the spectra of linear modes in e~e + plasma is meager. 
There exist only three modes Q| : purely transversal t mode and in general partially poten- 
tial It modes. Electric vector of t mode is perpendicular to both k and B . If k || B the 
dispersion relation of this mode is: 

Ut = kc(l-S), (5) 

where 5 = uij (Aug^ 3 )), and u p , ojb are plasma and cyclotron frequencies respectively 

(%ne 2 n p \ 1/2 e\B Q \ 

u p = , u B = , (6) 

\ m e J m e c 

angular brackets denote averaging over distribution function and (7) ~ j p . 

It mode has two branches lt\ and lt<i- Electric vectors of these modes are located in the 

plane formed by k and B . If k || B mode lt<i merges with t mode with spectrum (JSJ). The 

mode Iti propagating along Bo is longitudinal Langmuir wave (/). If u > kc their dispersion 

is HQ: 

^ = ^( 7 -) + 3^v(l-^), (7) 



and when uj « kc: 

oji » c [k - a(k - ho)] , (8) 

where k% = 2(7)^ /c 2 and a = (7)/2(7 3 ). 

From Eqs. (0) and (JHJ it follows that frequencies of I and t modes become equal at: 

^^^1 + ^<7 3 >J, (9) 

Dispersion relations of the waves propagating along B are presented in Fig. |21 

Mechanisms of wave generation in the magnetosphere e~e + plasma has been intensively 
studied (see, e.g., [ill)- The magnetosphere plasma can be unstable due to two reasons: a) 
existence of energetic beam; b) absence of symmetry and one dimensionality of the distri- 



bution function. Detailed analysis provides p] that the beam instability can develop only 
due to the drift of the particles in the curved magnetic field (so called " drift-cherenkov" 
resonance). Linear stage of this instability was studied in jl^ and nonlinear stage in J^. 
In this paper we will be interested in the second mechanism of instability - excitation at the 
anomalous Doppler-effect resonance. The condition of this resonance is jf|: 

uj - k\\v\\ - k ± u™ + — ~ 0. (10) 

7r 



Here hi and k± are projections of wave vector along and perpendicular to Bo, u™ = 
cv\\j r /uBRB is curvature drift velocity, Rb is curvature radius of magnetic field lines and 7 r 
is Lorentz factor of resonant particles. 

Detailed analysis provides that this condition can be satisfied only for energetic particles 



of the tale and beam 



21] , when the condition 2^5 > 1 is satisfied. Substituting Eq. (j2J) 



U) t » knC « — ^. (11) 



into (fTUj). for k± = we obtain: 



From this equation for typical pulsar parameters [22J, |23j it follows that the frequencies of 
generated t waves 

UJ t < cob- (12) 

It should be noted that transversal waves are strongly damped by cyclotron damping 
with the bulk particles. The condition of this resonance yields 

uj - k\\v\\ - k ± v% - — « 0. (13) 

This condition provides that cyclotron damping is effective for the waves with frequencies 

uj t > 2-i P u) B . (14) 

Linear mechanisms of the wave generation discussed in this section emerges transversal t 
waves. On its turn, nonlinear interactions of waves redistributes the energy among different 
modes and scales and as it will be shown bellow can lead to the formation of stationary 
spectra of wave turbulence. 

III. NONLINEAR PROCESSES 

For the development of weak turbulence theory the small parameter of weak turbulence 
is used. In the case of e~e + plasma this parameter is 

I F\ 2 

«1, (15) 



where E is electric vector of the wave. This condition implies that we consider the state 
when the energy density of plasma is much greater then the energy density of waves, except 
the parameter (JT5j) in the magnetosphere plasma there exist another two small parameters 
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that greatly simplify future analysis. The first small parameter is Up/u^ -C 1. The second 
one: 

„ (u — k\\v\\) 2 , . 

1p ,J V « 1 (16) 



allows us to use so-called drift approximation 

r~i 

As it is well known [25[ if turbulence is weak the dynamics is governed by three wave 
resonant processes. In e~e + plasma there exists additional limitations for realization of 
three wave resonant interactions. The absence of girotropy causes the fact that the terms 
proportional to the odd power of the charge do not contribute to the nonlinear currents. 
Another specific property of the magnetosphere plasma is that all the waves are propagating 
almost along the magnetic field. 



Breakup of I wave into two t waves was studied in |26j ]. and all possible three wave 
processes in j^). In this researches main attention was payed to possible breakup of I waves 
because it was believed that these waves should be the most unstable ones. Existence of the 
anomalous Doppler-effect resonance discussed in the previous Section leads to the conclusion 
that in fact the most interesting three wave processes are possible breakups of t waves that 
should be responsible for the formation of observed spectra of pulsars. 

For the analysis of nonlinear processes it is useful to introduce quantities and notations 
used in quantum physics. Usually plasmon occupation numbers N£ and wave amplitudes 



are introduced 



here 
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1 / d „ A \E, 



u k N x k = oo k < a x k a x k >= - —u 2 el ) i^-. (17) 
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e, A *e A 



£ 



^(kikj^ + uhij); (18) 



k o 

e A is polarization vector of \(lt,l,t) mode; Ej = ejE; is dielectric tensor; * denotes 
complex conjugation and angular brackets denote averaging over phases. 
Maxwell equations yield 

d t a k + iu^ = XajAa / dk 1 dk 2 du 1 du} 2 a k " 1 a\ 2 * 
x5(k - k x - k 2 )5(u -ux- uj 2 )Vx\x iX2 , 
where matrix element of the process A — > Ai + A2 



(19) 



x(^^ 2 £ A ) _1/2 Al (<9^W 2 £ A ) _1/2 A2 (T A |A 1 A 2 , 

LO=0J L . (J=CO. u 



(20) 



and (J\\x 1 x 2 is nonlinear conductivity tensor. 

If k± = 0, the process t — > t'+t" is forbidden due to the fact that the second order current 
j(2) ^ e 3 w^gn , ^ o, this process can take place but the matrix element of interaction 
is negligibly small |27(: Vt\t>t" ~ ^pA^i? ^ 

Another possible three wave process t — > t' + / was discussed in J^tJ without taking into 
consideration drift motion of the particles. In this case the frequency of generated Langmuir 
waves uji « k\c. If all the waves are propagating along Bo, matrix element of interaction is 

,2 



e ^ p loJtOJf 



VW«— ^ J — . (21) 

m e c V cj/ 

Consideration of the drift motion of particles makes this process even more intense. For 
simplicity let us consider the case when the waves are propagating along B . Let us assume 
that electric vectors of t and t' waves are directed along x and y axes respectively (see Fig. 
E} and B is parallel to z axis. E l x causes the drift motion of particles along y: 



< = &, (22) 

which on its turn generates nonlinear electric field: 

E = - c [u d , B'] . (23) 

It has to be noted that drift velocity is the same for electrons and positrons. Consequently 
the drift motion do not lead to current generation in the linear approximation. If the 
resonant conditions 

uj t = ui +u t ', k t = ki + k t >, (24) 

are fulfilled the beating of t and t' waves that generates longitudinal electric field (j2Bj) are 
in resonance with corresponding I wave. 

Analysis of the resonant conditions (}2"4*j) yields that I wave should have phase velocity a 
bit less then the speed of light. But this king of I waves in general can be strongly damped 
by collisionless Landau damping. In this case instead of three wave nonlinear interaction 
one has to consider scattering of t waves by plasma particles - so called wave-particle- wave 



interaction 



24, 



Analysis of standard pulsar parameters yield that in principle both cases are possible. 
In this paper we consider the case when / waves involved in resonance with t waves are not 
strongly damped by collisionless damping. 



Maxwell equations governing the dynamics in the case when all the waves are propagating 
along B yield: 

+4.^ = 0, (25) 



dt 2 dz 2 J " dt 
dE z . , SL 



()f +^Jz=0. (26) 

Nonlinear currents j^ L can be readily calculated as follows: substituting (J2~2~j) into Eq. 
()23|) and taking into account Faraday law we get: 

For we have 

H L = < Pe , (28) 

where perturbation of charge density p e can be determined from Poisson equation. 
Substituting in Eqs. (|23|) - (j2^j) electric field components as 

E z , ± (t) = E z Mt)e t(wt - kz) + c.c, (29) 

and taking into account Eqs. (|2T |) - (j2^j) as well as dispersions (0) and ^ we obtain for slowly 
varying amplitudes: 

1 WAV (3°) 



dE xkt 


i aoi ( k 1 ' 


dt 


= 2~B~o 


dE y k t , 


i cuji i 


f k l 


dt 


~2~Bo~ ' 




dE zkl 


i CLOl 


/ k f 


dt 


2 (l)B 





+ — J E zkl E* kt , (31) 



k 



(32) 



Assuming that there exist many waves with chaotic phases and introducing occupation 
numbers: 



E X ,yk\ 2 AT _ (Y)\E; 



N x , yk = N zk = w g /l ^ fc| , (33) 



and using standard technique [24, 25] we get for the matrix element of interaction: 



V t \w ~ -^-(coth' + u t ,k t )j ^r^\- ( 34 ) 

Comparison of Eqs. ()21|) and ()34|) yields that the drift motion of particles leads to important 
intensification of the considered process. 
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IV. STATIONARY SPECTRA 



In the case under consideration u)t,u>v 3> <^i- Consequently, the nonlinear interaction is 
local - only the waves with close frequencies can effectively interact. As it was discussed in 
Sec. 2 the generation of t waves takes place for uj <C ujb, whereas the cyclotron damping is 
effective for u > 2 r y p ujB- Nonlinear processes are responsible for the formation of stationary 
spectrum between these intervals of generation and silk. As it is known |lO| the stationary 
spectrum in the Kolmogorov interval is fully determined by the matrix element (JMj) an d 
linear dispersions of waves. There exist two stationary solutions: 

N kl = A 1 k' s ' d , N k2 = A 2 k~ s ~ d+1 '\ (35) 

where d is the dimension of the problem and s is the index of homogeneity of the matrix 
element. The first solution N k i = Aik~ s ~ d corresponds to constant flux of energy to the 
smaller scales, whereas N k2 = A 2 k~ s ~ d+1 ^ 2 corresponds to constant flux of plasmons to the 



larger scales 
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In our case d — 1 and 5 = 3/2. Therefore N k i = A\k~ b / 2 and N k2 = A 2 k~ 2 . Taking 
into account the linear character of t wave dispersion for the energy spectrum that is 
proportional to the observed intensity J w we have 

which of this spectra is realized in practice depends on the linear mechanisms of gener- 
ation and silk of the waves. In the case under consideration the waves with relatively low 
frequencies are excited at the anomalous doppler-effect resonance and waves with relatively 
high frequencies are damped at cyclotron resonance. In this case the spectrum E w \ ~ u~ 3 / 2 
is only possible. 

As it was mentioned above most of the pulsars has the spectral indexes — 1.5 > a > —2. 
Consequently, obtained theoretical result a = —1.5 seems to be in good accordance with 
the observations. 



V. SUMMARY 

Nonlinear three wave processes that can be responsible for the formation of spectra of 
radio pulsars is considered. Existence of strong magnetic field makes it necessary to take into 

9 



account electric drift motion of particles. In the framework of the weak turbulence theory- 
possible stationary spectra that can occur in Kolmogorov interval are calculated. Analysis of 
linear mechanisms of wave generation and silk allows to conclude that only possible spectral 
index is a = —1.5. Obtained result is in satisfactory accordance with the observational data. 
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FIG. 1: Particle distribution function / of magnetosphere plasma. Solid and dashed lines corre- 
spond to electrons and positrons respectively. 
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